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We consider the intermittent behavior of superfiuid turbulence in 4 He. Due to the similarity 
in the nonlinear structure of the two-fluid model of superfluidity and the Euler and Navier-Stokcs 
equations one expects the scaling exponents of the structure functions to be the same as in classical 
turbulence for temperatures close to the superfiuid transition T\ and also for T<T^. This is not 
the case when mutual friction becomes important. Using shell model simulations, we propose that 
for an intermediate regime of temperatures, such that the density of normal and superfiuid compo- 
nents are comparable to each other, there exists a range of scales in which the effective exponents 
indicate stronger intermittency. We offer a bridge relation between these effective and the classical 
scaling exponents. Since this effect occurs at accessible temperatures and Reynolds numbers, we 
propose that experiments should be conducted to further assess the validity and implications of this 
prediction. 



Introduction. Non-equilibrium systems are often char- 
acterized by relatively quiescent periods which are inter- 
rupted by rare but violent changes of the physical observ- 
ables. This phenomenon which is referred to as "inter- 
mittency" is very generic as it appears in many different 
stationary random processes. 

It is believed [1], that highly excited hydrodynamic 
turbulence possesses universal statistics in the inertial 
interval of scales L >• r >• r\. Here L is the of en- 
ergy input scale, e.g. due to instabilities of high velocity 
flows, and r\ is the viscous energy dissipation scale. The 
intermittent behavior of space homogeneous fully devel- 
oped turbulence of incompressible fluids can be studied 
in terms of the velocity structure functions [1] 

S p (r) = (\v(R,t)-v(R + r,t)\ p ) oc . (1) 

Here (...) denotes a time averaging. In the inertial inter- 
val, the functions S p (r) are scale invariant and they are 
characterized by the p-order scaling exponents £ p . 

Accepting Richardson's idea of a step-by-step energy 
cascade in the inertial interval and assuming that the 
energy flux over the scales e is the only relevant param- 
eter, Kolmogorov proposed from dimensional reasoning 
that 5p 41 (r) ~ (er) p/3 i.e. £ p = = p/3. In par- 
ticular, this means that S p (r)/[S2(r)] p ^ 2 should be r- 
independent, like in Gaussian statistics. However this 
scenario is not borne out by experiments [1]: the ob- 
served values £ oxp for p > 3 deviate down from p/3 and 
theratios5 p (r)/[5 2 (r)]P/ 2 increase with r. Therefore the 
probability of observing large values of small scale veloc- 
ity fluctuations (that mainly contribute to S p (r) with 
large p) significantly exceeds its normal K41 level, i.e. 
the small scale turbulent flow is strongly intermittent. 
In spite of fifty years efforts intermittency still requires 
deeper understanding. In particular there is no analytical 
framework to compute £ p (p). 

Turbulence in superfiuid 3 He and 4 He has emerged as 
a hot topic of interest merging the traditional fluid me- 
chanics with the low-temperature physics. Unlike classi- 



cal fluids, where the circulation around the vortices is a 
dynamical variable, in superfluids it is quantized to inte- 
ger values of k = 27rft/m. Here h is the Plank's constant 
and in is the mass of a superfiuid particle. This leads to 
the appearance of an additional length scale £, the mean 
inter-vortex distance. For motions with scales r » I 
superfluids can be described as two inter-penetrating en- 
tities with two different velocities: v s (inviscid superfiuid 
component) and v n (normal component), and tempera- 
ture dependent densities p s (T) and p n (T) = p — p s (T) 
interacting via the so-called mutual friction force [2]. 

One of the most interesting questions concerns the 
similarities and differences between superfiuid and clas- 
sical turbulence. The pioneering experiments of Mau- 
rcr and Tabeling suggest that intermittency would also 
be present in turbulent superfiuid 4 He [3]. This crucial 
observation implies that the existence of intermittency 
is not necessarily related to the exact structure of the 
Navier-Stokes equations and thus superfiuid turbulence 
offers an additional way to investigate its physical mech- 
anisms. Unfortunately, intermittency in superfluids has 
received far less attention than in classical turbulence [4] . 

The aim of this Letter is to compare intermittency in 
4 He-superfluid turbulence with turbulence in normal flu- 
ids. An important parameter characterizing superfiuid 
turbulence is the ratio p s /Pn- For T slightly below the 
phase transition temperature T\, when p s /p n <C 1, one 
can neglect the presence of the superfiuid component and 
the statistics of turbulent superfiuid 4 He is expected to 
be close to that of classical fluids. We also expect sim- 
ilar inertial range behavior of classical and superfiuid 
turbulence for T <C T\, when p n <C p s , due to the in- 
consequential role played by the normal component, see, 
e.g. [4]. Moreover (and less trivially) the intermittent 
scaling exponents £ p have to be the same in classical and 
low-temperature superfiuid turbulence. This is because 
the nonlinear structure of the equation for the superfiuid 
component is the same as the Euler equation, and dissi- 
pative mechanisms are believed to be irrelevant. 
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Based on the above reasoning, we focus in on the range 
of T ~ (0.9 0.8)T A , where /°s ~ Pn and compare it to 
the high and low T limits, when p s /p = 0.1 and 0.9, re- 
spectively. Our main tool is numerical simulations of su- 
pcrfluid turbulence in the two- fluid model [2] . We utilize 
a shell-model approximation, which is known to provide 
accurate scaling exponents in the classical case when the 
parameters are cleverly chosen (see, e.g. [5-7]). In par- 
ticular, shell models with the right parameters exhibit 
anomalous scaling exponents £p Um that agree well with 
experimental values ^ xp : £? um = 0.72 vs. £; xp = °- 70 ; 

£num = L26 vg . ^xp = 1 etc ^ 

The main result of the Letter is as follows: For tem- 
peratures such that p s /p < 0.1 and p s /p > 0.9, the scal- 
ing behavior of superfluid turbulence is practically the 
same as in classical fluids. However, for temperatures 
such that Ps/p ~ 0.5-^0.75, we discovered the emergence 
of a wide (up to 3 decades) interval of scales charac- 
terized by substantially different effective scaling expo- 
nents £ p ^ £p Um - For example, when p s = p n the value 
of £2 diminishes from its intermittent value ^J 1 ™ * 0-?2 

~ K41 

down to £2 ~ 0.67, close to the K41 value £ 2 = 2/3. At 
the same time, the exponents £ p for p > 3 deviate further 
from the K41 values: £ p < ^ um < i p . This fact can be 
considered as enhancement of intermittency in superfluid 
turbulence compared to the classical case. 

In addition, based on numerical observations that the 
probability distribution function (PDF) of the normal- 
ized shell energy is temperature independent, we suggest 
the relation between effective and classical exponents: 



P 



num 
2 



«) ~ 4 P 



0.026p, (2) 



which is in good agreement with observed values of £ p . 

Numerical procedure. Shell models of superfluid tur- 
bulence [8, 9] are simplified caricatures of the Navier- 
Stokcs and Euler equations in fc-representation: 



r d ,2 
dt + 
d 

dt + V » k 
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= nlk] + + r m , 

NL[u m ] = i{ak m+1 u m+2 u* m+l 



"t"^mUm+l u m -l 

a n , s n s (u s m - <J , 



-lU n 



^ km\ u m\ 2 



(3a) 
(3b) 
(3c) 

(3d) 



They mimic the statistical behavior of fe-Fourier com- 
ponents of the turbulent superfluid and normal velocity 
fields in the entire shell of wave vectors k m < k < k m+ i 
by complex shell velocity u*^ s . The shell wave num- 
bers are chosen as a geometric progression k m = fcoA" 1 , 
where m — 1,2, ... M are the shell indexes, and we have 



used the shell-spacing parameter A = 2 and ko = 1/16. 
In order to resolve a large range of scales we used a to- 
tal of M — 36 shells. This allows us to clearly resolve 
subregions with different scaling behavior. 

Similarly to the Navier Stokes equation, the NL[w m ] 
term in Eq. (3c) is quadratic in velocities, proportional 
to k and conserves (in the force-less, inviscid limit) the 
kinetic energy E = 5 X) m l u m| 2 , provided that a+b+c = 
0. We used the Sabra version [6] of NL[w m ] with the 
traditional (and physically motivated choice) b = c = 
—a/2, which ensures [5, 6] that £p Um are close to £p Xp - 

In Eq. (3a), v n is the kinematic viscosity of normal 
component. It is known that vortex reconncctions act 
as an energy sink in superfluids, which we accounted for 
(admittedly in a rough manner) by adding in Eq. (3b) 
an effective superfluid viscosity v s <C The details of 
this complicated and otherwise still not fully understood 
dissipation process go beyond the scope of this Letter. In 
simulations we used v s = 10~ 12 and v n = 10~ 10 . 

The coupling terms F^f [8, 9], given by Eq (3d), ac- 
count for the mutual friction between normal and super- 
fluid velocities. This dissipative force is proportional to 
parameters a s (T) — a and a n (T) — ap s / p n . 

The Gaussian (5-correlated in time forces f^f were ap- 
plied only for to = 1,2. Eqs (3) were solved using 4-th 
order Runge-Kutta method. The chosen values of p s /p 
and corresponding values of T/T\ and a (interpolated 
using the experimental data in [10]) are shown below: 
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Results and discussions. First of all, we study the 
scaling behavior of the velocity structure functions (1), 
which in the shell models can be written as 

S s p ' n (k m ) = (fcT) , C p {k m ) = kk S p (k m ). (4) 

We defined for convenience the structure functions 
C p (k m ) compensated by the classical anomalous scaling. 

Figure 1A shows C-2(k m ). The almost horizontal lines 
for p s /p = 0.9 and 0.1, indicate the usual intermittent 
scaling behavior. However for p s /p — 0.5 one detects 
the emergence of a wide enhancement interval with the 
energy accumulations and where the scaling exponent is 
close to £2 — 2/3 [11]. This is clearly seen in the inset, 
where ^(fcm) is compensated with K41 scaling. At first 
glance, the scaling with £ 2 — £2 " =2/3 can be consid- 
ered as the normal K41 scaling in superfluids. In this 
case one would expect that the PDF of normalized en- 
ergy m to- shell P n > s [|<f | 2 /S£' s (fc TO )] should collapse for 
all to. This is, however, NOT the case here. Figure IB 
shows that the PDFs for different to are in fact very dif- 
ferent from each other. The much longer tails associated 
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FIG. 1: Panel A: The compensated structure functions C2 ,s (fc m ) for p s /p = 0.1,0.5 and 0.9, normalized by its mean value over 
shells with m = 7 — 14. Inset: fc 2//3 5*2 (k m ) normalized by its mean value over m = 15 — 20 for p s /p = 0.5. Solid lines - normal 
component, dashed lines - superfluid components. The upper axis is marked with the shell numbers, the lower axis shows the 
corresponding wavenumbers. Panel B: Demonstration of the temperature-independence of the PDF's for shells m = 14 and 
m = 20. Panel C: The normalized compensated C2' s (k m ). Inset: The collapse of k s S4,(k)/[S'2(k)] 2 , compensated using classical 
exponents S = 2£ 2 mm — ^" um ~ 0.184. The horizontal black dashed lines serve to guide the eye only. Color code and line labeling 
is the same in all panels. Vertical dot-dashed lines in Panels A and C denote the intervortex scale l/£ (see [8]). 



with larger values of k m reveal strong intermittency ef- 
fects. This is clearly seen in Fig. 1C where C^'^fcm.), 
is almost flat for p s /p — 0.9 and 0.1, but show a dif- 
ferent scaling for p s /p — 0.5 with an effective exponent 

— K41 

£4 « 1.21 < £J um = 1.256 < £ = 4/3. The next or- 
der structure functions, Se{k), S$(k),.. demonstrate the 
same type of behavior where the effective exponents £ p 
deviate even further away from the K41 values than the 
intermittency exponents £p Um in classical fluids. 

Another important observation in Fig. IB is that the 
PDFs are practically temperature independent. Then 
knowing £ p um and £2, we can predict with very good 

accuracy all the other effective exponents £ p . Indeed, 
comparing the following expressions for Sp ,s (k m ) 

/>oo 

S£ s (k m )= / V^[\u^\ 2 /Sr(k m )}\u^rdu^ s 
Jo 

in the regions of different scaling, one arrives after some 
simple algebra to Eq. (2). An immediate consequence of 
Eq. (2) is that ratio 5 , p (fc)/[5 , 2 (A:)] p / 2 is the same in the 
regions with £ p um and £ p scaling. This prediction fully 
agrees with observations (see e.g. inset in Fig. 1C for 
p = 4) and thus lends strong support to Eq. (2). 

To quantify the intermittency enhancement we dis- 
play in Fig. 2 the T- and p s /p-dependence of £2 to- 
gether with the dimensionless magnitude $2 of this effect, 
defined as ^2 = log 2 {[maxC2(A: m )/(C2(fc m ))]}. Here 

m 

(C2(k m )) is the mean value of C'2(k m ) on the plateau. 
Denote by m± the left and right edges of the enhance- 
ment interval and its extent as k m+ /k m _ — 2 Am , Am = 
(m + — m_). Then, using data in Fig. 2 and comput- 
ing the ratios 4 r 2 (T)/[£|- um - ^(T)] = Am for all T, we 
found that Am ~ 8 -7- 10 practically independent of T. 
In other words the extent of the enhancement interval 
k m+ /k m - = 2 Am ^ 250 1000. 

As expected, significant intermittency enhancement 



takes place for p s ~ p n [or more precisely p s ~ (0.4 
0.8)/? - see Fig. 2] when the two-fluid model (3) differs 
mostly from a single classical fluid. Normal and super- 
fluid components in Eqs. (3) are coupled by the mu- 
tual friction force F m oc ot(u s m — it£j which causes the 
energy exchange between normal and superfluid com- 
ponents and the energy dissipation. Thus these two 
phenomena are the only possible reasons for the ob- 
served intermittency enhancement and both disappear 
in case of full velocity locking u s m {t) = Veloc- 
ity locking can be statistically quantified by the cross- 
correlation coefficient of normal and superfluid veloci- 
ties: K m = + <,<»*])/([«£<. + u m<T])- For 
fc(k m ) = 0, the superfluid and normal components are 
statistically independent subsystems, while for /C(fc m ) = 
1 they are completely locked to each other. Clearly only 
for < JC(k rn ) < 1 we can expect intermittency en- 
hancement. A result of [13] for the correlation coeffi- 
cients can be reformulated for shell models as: K, m ~ 

2[a n ^, m + « g ^,J/[5I, m +5 2 n ,J 
a s + a s + k m [v n k m + 5| m +S 2 \ m ] / VTW) 
Here SJm = (k m ). Analysis of this equation shows 
that the and are locked when k m i] <C 1 and be- 

T/T, 
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FIG. 2: T/T\- and p 3 (T)/p-dependencies (upper and lower 
axes) of the enhancement magnitude ^2, left axis, and the 
exponent £2, right axis. The maximum of ^2 ~ 0.7 and min- 
imum of £ 2 « 0.65 are reached at T w 0.865T\ (p s « 0.6p). 
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come uncorrelated when k m ~ (ap/p n ) /?7. As one sees 
in the Table, for all T in our simulation ap/p n ~ 0.5H-1.2. 
Therefore one should expect velocity locking for all T and 
for all k m r\ < 1 including the enhancement range. This 
prediction agrees with our results. 

This means that the origin for the variations of the 
scaling is not in the enhancement range but rather in 
the vicinity of the viscous cutoff, where k + <~ k v . This 
conclusion is strongly supported by our simulations with 
different T, M = 28 + 46 and v n /v s ~ 10 2 -j- 10 6 , in all 
of which the extent of the enhancement interval remains 
unchanged, Am ~ 8 10, and always k + ^ k v . 

One concludes that a key role should be played by the 
energy dissipation due to the mutual friction D m — 
ap s n s \uj n — u"J 2 , as follows from Eqs. (3). Clearly, 
D™ F — > for T —¥ T\ (because p s -> 0) and for T -> 
(because a — >• 0) and has a maximum at T ~ 0.80-^0.85, 
i.e. in the range, where the intermittency enhancement 
is observed. Moreover, in this range, the total mutual 
friction dissipation ^2 m D m is close to (or even larger 
than) the total viscous dissipation ^nJ2m^m^2,m- im ~ 

MF 

portantly, D m has a maximum at k m , close to that of the 
viscous dissipation. As a result, for T ~ (0.80 -j- 0.85)T\, 
the mutual friction between supcrfiuid and normal com- 
ponents markedly drives the system away from the simple 
one-fluid behavior primarily around k ~ k v , leading to a 
sharper decay of the energy spectrum S^'^,. I n turn, this 
may give rise to a bottleneck effect leading to an energy 
accumulation near k v . This phenomenon shares similari- 
ties with the bottleneck in classical fluids [14] but differs 
in important aspects: e.g. as in classical fluids the width 
of the energy accumulation range is almost independent 
of underlying mechanism and only its magnitude varies. 
In our case, however, the extent of the enhancement in- 
terval is much larger, about 3 decades, compared to one 
decade in classical fluids. A possible reason for such 
an extended enhancement interval is the energy transfer 
from supcrfiuid to normal component, caused by mutual 
friction, and it is also notable in the same T range. 

A comprehensive study into the quantitative aspects 
of this enhancement of intermittency is beyond the scope 
of the Letter. It will require an extensive analysis of the 
time dependencies of u s m {t) and in order to clarify 

the role of the mutual friction on the statistics of outliers; 
the main contributors to the fat tails of the PDFs which 
are responsible for the intermittency [15]. 
Summary. Our shell-model simulations of the large- 
scale superfluid turbulence in 4 He indicate that for high- 
and low- T- limits, when p s < O.lp and p s > 0.9p anoma- 
lous scaling exponents of superfluid turbulence are prac- 
tically the same as their classical counterparts, as fol- 
lows from our theoretical reasoning and agrees with nu- 
merical and laboratory observations [4] . The shell- model 
simulations have allowed us to uncover considerable T- 
dependent variations of the apparent scaling over the 
enhancement interval of scales about 3 decades when 



Ps ~ Pn- This scaling can be characterized by the ef- 
fective scaling exponents £ P (T) that deviate further from 
the normal K41 values as compared to the exponents of 
classical fluids. Furthermore if the entire inertial inter- 
val is shorter than the enhancement interval, only the 
effective scaling is observed. We also suggested Eq. (2) 
relating £ P (T) to £ p . 

While our interpretation indicates the existence of a 
form of the energy bottleneck in the two-fluid system, 
the phenomenon definitely calls for further investigation. 
DNS of supcrfiuid turbulence oriented towards the study 
of the predicted T-dependent scaling over an interval of 
about three decades is nowadays realistic and appears 
as particularly timely. Even more importantly, the tem- 
perature range T ~ (1.7 -j- 1.9) K, where we forecast a 
change in the scaling behavior is accessible to laboratory 
experiments. We hope that the results presented in this 
Letter will inspire new numerical and laboratory studies 
of intermittency in superfluid turbulence. 
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